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. Abstract. We present recent multi-object spectroscopy of globular clus- 

O ' ters in the Virgo gEs M49 and M87. In M49, we have a total of 144 con- 

, firmed clusters out to 8 arcmin radius (~ 6 Re// or 35 kpc). We find that 

P3 ' the blue (metal-poor) clusters have both a higher velocity dispersion and 

^ . rotation than the red (metal-poor) clusters. For the metal-rich population 

we place an upper limit of {y/(T)proj < 0.34 at 99% confidence. We cal- 
culate the velocity dispersion as a function of radius, and show that this 
is consistent with isotropic cluster orbits and the M49 mass distribution 
^ ■ determined from X-ray data. For M87, we combine new CFHT data with 

00 I previous data to obtain a total sample of 278 clusters out to 10 arcmin 

■ radius {'^ 45 kpc). We find a similar global rotation for the metal-poor 

and metal-rich clusters of 160—170 km/sec. Beyond ~ 2 Re// (15 kpc), 
^ ■ both the metal-poor and metal-rich clusters appear to rotate about the 

photometric minor axis. The combined cluster sample is consistent with 
O ■ isotropic orbits, but when considered separately, the metal-poor clusters 

show significant tangential bias of f3ci — —0.4, while the metal-rich clus- 
CL|' ters show a radial bias with Pel — -1-0.4. In both galaxies, the metal-rich 

Q , and metal-poor clusters share different kinematics, but there is no clear 

i-{ ' preference for any one formation scenario. 



X 



1. Introduction 

There is still considerable controversy about how and when elliptical galaxies 
form. As survivors from the earliest epochs of galaxy formation, globular clusters 
(GCs) are incredibly useful probes of their host galaxies. Of great recent interest 
has been the finding of himodality in the GC colour distributions in many ellip- 
tical galaxies. Several models have been put forward to explain this bimodality, 
including spiral-spiral mergers (Schweizer 1987; Ashman & Zepf 1992), multi- 
phase collapse (Forbes, Brodie, & Grillmair 1997), and hierarchical accretion 
(Cote, Marzke, & West 1998). 

The kinematics and ages of metal-poor and metal-rich GCs around ellip- 
ticals will be key to distinguishing between these different scenarios, hence the 
importance of cluster spectroscopy. However, this spectroscopy is difficult, given 
the faintness of GCs in even the nearest ellipticals. It has only been with the 
advent of multi-object spectrographs on 4m and 8— 10m telescopes that real 
progress has been made in this area. However, the sample of ellipticals with 
signficant numbers of GC spectra is still very small, limited to giant elliptical 
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(gE) galaxies such as Centaurus A, NGC 1399, and M49 and M87 in Virgo. Here 
we present recent multi-slit spectroscopy of GCs in the latter two galaxies. 

2. M49 Globular Clusters 

With K. Ashman, M. Beasley, D. Geisler, D. Hanes, R. Sharpies, and S. Zepf. 
See Sharpies et al. (1998) and Zepf et al. (2000) for further details of the 
M49 GC kinematics, and Beasley et al. (2000) regarding the GC ages and 
metallicities. 

2.1. Sample Selection and Observations 

Object selection was based on Washington photometry of M49 chistcr candi- 
dates from Geisler, Lee, & Kim (1996). To improve the GC yield, we imposed 
magnitude and colour cuts: 19.5 < V < 22.5, 0.5 < C— Ti < 2.2. Cluster spectra 
were obtained at the WHT with LOSS in 1994, and at the CFHT with MOS in 
1998. Spectra covered the region from 3800—6000 A with a resolution of 3—6 
A, giving velocities good to 50—100 km/sec; exposure times were 3—4 hours 
per mask. The combined cluster sample from the two runs (also including some 
GCs from Mould et al. 1990) is 144; of these, 93 are blue/metal-poor and 51 are 
red/metal-rich (dividing at C— Ti = 1.625, or [Fc/H] ~ —1.6). We suffered only 
15% contamination from non-clusters, and our overall completness was 65%. 
Our sample extends out in radius to ~ 8 arcmin (~ 6 Re//, or ~ 35 kpc for 
D=15 Mpc). Cluster velocities were obtained from cross-correlation with radial 
velocity stars. 

2.2. Results 

Kinematics Figure 1(a) shows the GC rotation and velocity dispersion radial 
profiles, for the total sample, and for the blue and red GCs separately. These 
profiles were obtained by smoothing with a Gaussian kernel with a = 100" . 
The blue GCs have a significantly larger velocity dispersion than the red GCs 
(300—350 km/sec compared to ~ 200 km/sec), as already found by Sharpies 
et al. (1998). Figure 1(a) shows that the M49 GCs have a slowly declining 
velocity dispersion profile, seen in both the full and blue/red samples. There is 
little evidence for rotation in the red GCs (but some evidence that the rotation 
increases outwards), while the rotation is roughly constant for the blue GCs. At 
all radii, the rotation is much smaller than the dispersion, and for the red GCs 
we can set an upper limit of {y/(T)proj < 0.34 at 99% confidence; thus, rotation 
is not dynamically important for the M49 GCs. 

M/L Ratio and Orbital Anisotropy We use the velocity dispersion and lumi- 
nosity density profiles, together with the Jeans equation, to estimate the M49 
mass distribution; we assume isotropic orbits and spherical symmetry. Figure 
1(b) shows the M49 mass vs. radius; the lines represent our cluster data, and the 
points with errorbars are the masses inferred from ROSAT X-ray data (Irwin & 
Sarazin 1996). There is reasonable agreement between the two mass determi- 
nations, suggesting that both are roughly correct, and that our assumption of 
isotropic cluster orbits isn't too far wrong. We find a total M49 mass in excess 
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Figure 1. (a) (left): Smoothed rotation and velocity dispersion pro- 
files for the M49 GCs. Top panels: combined sample; Middle panels: 
blue (metal-poor) GCs; Bottom panels: red (metal-rich) GCs. Dotted 
lines show the la bootstrapped uncertainties. 

(b) (right): M49 mass vs. radius. Solid line: mass determined from 

the cluster radial velocities, assuming isotropic orbits; Dotted lines: 
la bootstrapped uncertainties; Points with error bars: masses from 
ROSAT X-ray data (Irwin & Sarazin 1996). 

of 10^^ M0, and a M/L ratio at least five times greater at 30 kpc than at a few 
kpc, confirming the existence of a substantial dark matter halo in M49. 

3. M87 Globular Clusters 

with P. Cote, D. Ceisler, D. Hanes, G. Harris, J. Hesser, D. McLaughlin, and 
D. Merritt. See Hanes et al. ( 2001 ) and Cote et al. ( 2001 ) for further details. 

3.1. Observations and Final Sample 

Object selection was based on Washington C, Ti imaging (Gcislcr, Lee, & Kim 
2001), with 0.8 < C-Ti < 2.35. Spectra were obtained at the CFHT in 1996 
with MOS, with a dispersion of 3.6 A/pixel. We used a filter centered at 5100 A 
with a bandwidth of 1200 A to allow us to place up to 100 slits per mask with- 
out spectral overlap. Galactic GCs were observed as radial velocity templates. 
We obtained velocities for 109 M87 GCs, 54 of these new measurements. By 
combining with previous work (most importantly, 205 velocities from Cohen k, 
Ryzhov 1997), we obtain a final sample of 278 bona-fide M87 GC velocities. Of 
these, 161 are classified as "blue" (C— Ti < 1.4) or metal-poor ([Fe/H] < —1.15) 
and 117 as "red" or metal-rich (C-Ti > 1.4, [Fe/H] > -1.15). 

3.2. Results 

Cluster Kinematics We investigate the global rotation and velocity dispersion 
of the M87 GCS by fitting our velocities with the function Vp{0) = Vgys + 
(Qi?)sin(G — Go), where Qq is the GCS rotation axis, and fiR is the rotation 
amplitude. When we average over all radii, we find similar rotation, dispersion, 
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and position angles for the full sample, and for the blue and red clusters sepa- 
rately: ^IR ~ 160—170 km/sec, a dispersion of 360—400 km/sec, and a position 
angle between 60—75°. 

When we look at the kinematics as a function of radius, however, it starts to 
get interesting. Figure 2(a) shows the cluster rotation amplitude vs. radius for 
the total sample, and for the red and blue clusters separately, while Figure 2(b) 
shows the cluster best-fit position angle vs. radius. These (smoothed) curves 
were obtained by sliding a bin of radial width AR= 90" (~ 6.5 kpc) through 
the cluster datapoints, and solving for the kinematics in each bin. From Figure 
2(a), we see that the rotation is roughly constant with radius, but there is some 
evidence that the rotation increases at larger radius, especially for the blue GCs. 
Figure 2(b) shows that, at large radius, the blue GCs rotate conventionally 
around the minor axis of M87, but for R < 16—18 kpc, they appear to rotate 
about the major axis! Similar behavior is seen in the red GCs, though not as 
pronounced. This flip in the rotation axis may be evidence for a past major 
merger. 





Figure 2. (a) (left): Projected rotation velocity vs galactocentric 

radius for the M87 GCS (filled circles), with 68% and 90% confidence 
bands (thin and thick solid lines). The smooth, dotted curves repre- 
sent the circular velocity of the M87/Virgo potential from the model 
in McLaughlin 1999. Upper panel: combined sample; Middle panel: 
metal-poor GCs; Bottom panel: metal-rich GCs. 

(b) (right): Same as Figure 2(a), but now plotting the projected az- 
imuth of the rotation axis of the M87 GCS vs galactocentric radius. 
The lower and upper horizontal lines show the position angles of the 
M87 major and minor axes. 



Dynamical Models and Cluster Orbits Our dynamical approach is to assume 
a total M87 mass profile, as given by McLaughlin (1999). We plug this mass 
profile, the deprojected globular cluster surface density profile, and the velocity 
anisotropy I3ci{r) = 1 — a^/a^ into the Jeans equation and solve for the radial 
velocity dispersion o"r(r). A deprojection then gives the projected velocity dis- 
persion profile o"p(r), which can then be compared with the observed profiles. 
Thus, we input various values of/? [specifically (3 = —99 (strong tangential bias). 
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—0.4 (moderate tangential bias), (isotropic), +0.4 (moderate radial bias), and 
0.99 (strong radial bias)], and see which value fits the data best. 

Figures 3(a-c) show the results of this, for the combined sample, and for the 
metal-poor and metal-rich GCs separately. From Figure 3(a), the M87 GCS as 
a whole has an almost perfectly isotropic velocity ellipsoid. However, the metal- 
poor GCs appear to have a modest tangential bias, with P^l ~ —0.4 at small 
radii, while the mctal-rich GCs appear to have a slight radial bias of roughly the 
same magnitude, ~ +0.4. Thus, the metal-poor and metal-rich GCs appear 
to have quite different orbits. 




R/kpc 



R/kpc 



R/kpc 



Figure 3. (a) (left): Upper Panel: Velocity dispersion profile for 
the combined sample of M87 GCs. Points: smoothed data (the thin 
solid curves are 90% confidence intervals); Thick solid line: predicted 
velocity dispersion profiles for isotropic orbits (/3= 0); Dashed lines: 
/3=0.99, 0.4, —0.4, and —99, from top to bottom. Lower Panel: similar 
to upper panel, except now for the aperture velocity dispersion profile. 

(b) (middle): Same as Fig 3(a), except for the metal-poor GCs. 

(c) (right): Same as Fig 3(a), except for the metal-rich GCs. 



4. Discussion: Comparison of M49 and M87 

What can we learn by comparing the M49 and M87 GCSs? 

• In M49, there is little or no rotation in the red GCs, and ^ 100 km/sec 
rotation for the blue GCs. In M87, there is a similar rotation of ~ 170 
km/sec for both red and blue GCs. 

• In M49, there is little or no evidence that the rotation or velocity disper- 
sion is rising with radius, whereas in M87 there is a suggestion that both 
the rotation and dispersion are increasing outwards. This is most likely 
because M87 is at the center of Virgo, and the M87 GCs are tracing the 
transition to the Virgo cluster potential. See Cote et al. (2001) for more 
discussion on this point. 

• In M49, V/o- ~ 0.3/0.1 for the blue/red GCs, while in M87 Y/a ~ 0.45 
for both subpopulations. 
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• The combined samples of GCs in M49 and M87 arc both consistent with 
isotropy. In M87, the blue GCs appear to have a tangential bias, while the 
red GCs have a radial bias. 

Thus, in both galaxies, the metal-poor and metal-rich GCs have different 
kinematics (rotation, velocity dispersion, and orbital anisotropy), but these dif- 
ferences are not consistent between the two galaxies. It is difficult to think of 
a single formation model that can explain the properties (kinematic and other) 
of both GCSs. There is obviously a huge need for further kincmatical samples 
in a wider range of galaxies. Luckily, we will soon have fabulous new multi-slit 
spectrographs on 8m telescopes: CMOS on Gemini and VIMOS on the VLT. 
The next few years will be very exciting. 
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